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O ■ ABSTRACT 

?^ ! We discuss the effects of resonant spin-ffavor precession (RSFP) of Majo- 

\ rana neutrinos on lieavy element nucleosynthesis in neutrino- heated supernova 

^ \ ejecta and the dynamics of supernovae. In assessing the effects of RSFP, we 

■ exphcitly include matter-enhanced (MSW) resonant neutrino flavor conversion 

effects where appropriate. We point out that for plausible ranges of neutrino 
magnetic moments and proto-neutron star magnetic fields, spin-flavor conversion 
of v-r (or z/^) with a cosmologically significant mass (1-100 eV) into a light 
could lead to an enhanced neutron excess in neutrino-heated supernova ejecta. 
This could be beneficial for models of r-process nucleosynthesis associated with 
late-time neutrino-heated ejecta from supernovae. Similar spin-flavor conversion 
of neutrinos at earlier epochs could lead to an increased shock reheating rate and, 
concomitantly, a larger supernova explosion energy. We show, however, that such 
increased neutrino heating likely will be accompanied by an enhanced neutron 
excess which could exacerbate the problem of the overproduction of the neutron 
number = 50 nuclei in the supernova ejecta from this stage. In all of these 
scenarios, the average z/g energy will be increased over those predicted by super- 
nova models with no neutrino mixings. This may allow the SN1987A data to 
constrain RSFP-based schemes. 
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I. Introduction 



In this paper we examine the effects of Resonant Spin-Flavor Precession 
(RSFP) of Majorana neutrinos on heavy element nucleosynthesis in neutrino- 
heated supernova ejecta and the dynamics of supernovae. Massive neutrinos that 
possess a transition magnetic moment could experience a resonant conversion 
of spin and flavor in the presence of magnetic fields and matter This 
RSFP effect is similar to the well-known Mikheyev-Smirnov-Wolfenstein (MSW) 
mechanism 0. In fact, in the environment above the neutrinosphere in post-core- 
bounce supernovae, both the MSW and RSFP processes may operate. Reference 
discusses the effect of RSFP-induced neutrino flavor conversion on supernova 
shock reheating. Here we extend their discussion by considering the effects of 
RSFP on nucleosynthesis in late-time neutrino-heated supernova ejecta and gen- 
erally including the effects of both RSFP and MSW conversion. 

Supernovae have long been considered as a promising site for heavy element 
nucleosynthesis . About one half of the elements with mass number A > 70 in 
nature are believed to be made by the rapid neutron capture process, or r-process 
for short. In the r-process, neutron captures occur much faster than typical 
/3-decays. Recent calculations suggest that the high-temperature, high-entropy 
region which would be formed above the proto-neutron star several seconds after 
the bounce of the core in a type-II supernova is a promising site for r-process 
nucleosynthesis [§], 0, 

This site is sometimes referred to as the "hot bubble," since material there is 
heated by absorption of neutrinos and antineutrinos which are emitted from the 
neutrinosphere near the surface of the hot proto-neutron star. Close to the neu- 
trinosphere the temperature is high enough that all strong and electromagnetic 
nuclear interactions are in equilibrium (nuclear statistical equilibrium, or NSE). 
As the material above the neutrinosphere expands due to neutrino heating, its 
temperature and density decrease. When the temperature drops below about 
0.5 MeV, the material outflow rate (or expansion rate) becomes faster than the 
rates for nuclear reactions. At this point, the material freezes out of NSE. As 
the material further expands above this nuclear freeze-out point, rapid neutron 
captures onto the existing seed nuclei may occur after an alpha-rich freeze-out of 
the charged-particle reactions. 

In order for the r-process to occur, the material certainly has to be neutron 
rich at the nuclear freeze-out position. The neutron-to-proton ratio above the 
neutrinosphere is determined by the following reactions, 

Ue + n — > p + e~, (la) 
i^e+P n + e"*". (lb) 

Because material in the surface layers of the proto-neutron star consists mostly of 
neutrons, have a larger opacity than z/g [cf. Eqs. (^) and (^)]. Consequently, 
Ue decouple deeper inside the proto-neutron star where it is hotter, and corre- 
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spondingly have a higher average energy than z/g. In turn, this average energy 
hierarchy favors the rate of the process in Eq. over that in Eq. ([Ta|) . These 
arguments suggest that neutron-rich conditions which are conducive to r-process 
nucleosynthesis will obtain at the nuclear freeze-out point and beyond, in the 
region where the r-process may occur. 

Reference discusses these weak reaction issues and the connection between 
the flavor mixing of neutrinos with cosmologically significant masses and the 
conditions necessary for heavy element nucleosynthesis in supernovae. They show 
that resonant flavor conversion of a or z/^ with a mass of 1-100 eV into a 
light Ve may preclude heavy element nucleosynthesis in the hot bubble unless 
the vacuum mixing angle satisfies sin^2^<10~^. Because and and their 
antineutrinos lack the charged- current reactions similar to those in Eqs. (0) and 
(JTBp, they decouple deepest inside the proto- neutron star and have the highest 
average energy. Consequently, a significant amount of z/^(^) Ve transformation 
results in z/g with an average energy higher than that of z/g. This enhances the 
rate of the process in Eq. ( p^ and drives the neutrino-heated supernova ejecta 
proton-rich. 

If the mass of z/^ (or z/^) is between about 1 and 100 eV, then matter-enhanced 
resonant MSW fiavor conversion could occur in the region between the neutri- 
nosphere and the radius where the weak reactions in Eqs. (p!aD and ( p!B|) freeze 
out (very near the nuclear freeze-out position in most supernova models). If 
average energy or z/^ are converted at resonance with greater than about (25- 
30)% efficiency, then the hot bubble will be driven proton- rich and r-process 
nucleosynthesis at this site will be impossible. 

With no such fiavor conversion, the neutron excess required for r-process nu- 
cleosynthesis may or may not be obtained currently in models of neutrino-heated 
supernova ejecta. Although the conditions of entropy, electron fraction [a mea- 
sure of the neutron-to-proton ratio and the neutron excess, cf. Eq. (0)], and 
expansion rate as computed in some hydrodynamic models [§] provide the requi- 
site neutron-to-seed ratio for the r-process to occur, these models have left out 
key physics input that can wreck the r-process [|10|. Besides, simple wind model 
arguments suggest that the neutron-to-seed ratio obtained in these hydrodynamic 
models is unrealistically large on account of their high entropies at late times. 
Wind models suggest an entropy roughly half of that obtained in some hydrody- 
namic models, and this would imply a neutron-to-seed ratio too low to allow the 



production of the heaviest r-process nuclides [11, 12 . 



Any effect that could /ott;er the electron fraction Y,, (i.e., raise the neutron-to- 
seed ratio) in these models would be most welcome. In fact, it is even conceivable 
that r-process nucleosynthesis in late-time neutrino-heated supernova ejecta will 
be impossible unless there is some new physics which has the effect of raising the 
neutron excess [0, |T^. However, any effect that significantly lowers Y^. at late 



times must not also do so at the early times characteristic of shock reheating, 
lest the neutron number = 50 nuclei be overproduced |T3|, |T^. This shock 
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reheating epoch occurs roughly ~ 0.1-0.6 s after core bounce, as opposed to the 
epoch where the r-process might take place in neutrino-heated ejecta at ~ 3-20 
s after core bounce. 

In what follows we show that (or z/^) with vacuum masses in the range 1-100 
eV could convert in principle into Ve by RSFP above the neutrinosphere. In turn, 
this could lead to the enhancement of the neutron excess in the "hot bubble." In 
order to have the conversion i^Tdi) ^ neutrinos must be of Majorana type and 
have a finite transition magnetic moment. The operation of RSFP also requires 
a (probably large) magnetic field around the proto-neutron star. If RSFP did 
occur, Ve would become much more energetic than z/g and this would enlarge 
the neutron-to-proton ratio as outlined above. This implies that RSFP could be 
beneficial to the production of the heavier r-process elements in supernovae, in 
contrast to the case of resonant MSW flavor conversion alone as discussed in Ref. 
0. For significant RSFP effects to occur, the required product of the neutrino 
transition magnetic moment and magnetic field B near the neutron star will 
be shown to be of order 1 in units of Bohr magneton times gauss (/i^ ■ G). 

Before entering into a detailed discussion of RSFP in supernovae, let us de- 
scribe the present upper limits on the neutrino magnetic moment from laboratory 
experiments and astrophysical arguments. The present upper bound on the neu- 



trino magnetic moment from v^e scattering experiments is |T7 

/i^ < 1.9 X IO^^Vb- (2) 

This bound applies to the direct or transition magnetic moment of Dirac neutri- 
nos, as well as the transition magnetic moment of Majorana neutrinos. 

A stronger limit for the neutrino magnetic moment can be derived from well- 
known arguments against excessive cooling of red giant stars. A finite neutrino 
magnetic moment would enhance the plasmon decay into vv pairs in the stellar 
interior, resulting in excessive cooling. The most severe constraint is derived by 
estimating the critical mass for a helium flash in red giant stars. This gives the 
bound 



/i, < 3 X 10-^Vb (3) 

for the transition magnetic moment of Majorana neutrinos. 

A magnetic moment of order (10~^^-10~^°)/is is very large in the context of 
the small neutrino masses we might consider here. In general, this is because 
the diagram that generates a magnetic moment of order 10~^^/xs or larger with 
the photon line removed also induces a large neutrino mass. There are many 
(successful but not compelling) attempts at constructing a mechanism to induce 
a large neutrino magnetic moment, while keeping the masses of neutrinos small 



r9| . Since these issues are very speculative, we will assume here that a neutrino 
transition magnetic moment of order 10~^^/xb is plausible. However, we will 
see that the crucial quantity governing the effects of RSFP in supernovae is the 
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product of the neutrino transition magnetic moment and the magnetic field. 



II. Resonant Neutrino Spin-Flavor Precession in Supernovae 

In this section we discuss some general features of RSFP in supernovae. We 
also examine the case where both matter-enhanced MSW fiavor conversion and 
RSFP occur in the region above the hot proto-neutron star. Both Dirac and 
Majorana neutrinos can have RSFP, so long as the transition magnetic moment 
exists. The mechanism of RSFP is essentially the same for both cases. However, 
the implications of RSFP for supernova heavy element nucleosynthesis and/or 
explosion dynamics will be very different for the two cases. This is because RSFP 
for Majorana neutrinos occurs between two active neutrinos, whereas RSFP for 
Dirac neutrinos occurs between active and sterile neutrino states. In this paper 
we shall consider only the Majorana neutrino case. 

The Lagrangian which describes the magnetic moment-mediated interaction 
between Majorana neutrinos and the electromagnetic field, Fap, is given by 



where {fi)ab is the magnetic moment matrix with a, b = e, fi, t or 1, 2, 3 for the 
flavor or mass eigenstate bases, aa/3 = (i/2)[7a,7^] with 7^ the Dirac matrices, 
and C denotes the operation of charge conjugation. From the requirement for 
CPT invariance, the diagonal elements of the magnetic moment fiaa vanish, and 
consequently a transition magnetic moment is the only possibility for Majorana 
neutrinos. In the case of a finite neutrino transition magnetic moment, the pres- 
ence of magnetic fields can facilitate the transformation UaL {^bL}'" 7^ b) or 
vice versa. Since {ybi)^ is generally termed Vh (antineutrino state for z/fe) and is 
right handed, we can describe the Va ^ i^b {O' ^ b) transformation as a "spin- 
fiavor" precession (or conversion). Except for the interaction in Eq. (^, we will 
assume here that neutrinos possess only standard electroweak interactions with 
matter. 

Among the several conceivable channels of spin-fiavor precession, Vrip.) ^ is 
the example we will consider in what follows. Motivated by the average neutrino 
energy hierarchy discussed in the last section, we anticipate that this channel of 
RSFP may give the most significant effect on the electron fraction. We will show 
that the spin-flavor conversion of (or z/^) with masses in the range 1-100 eV 
into a light v^. can result in important effects on the parameters which determine 
heavy element nucleosynthesis and/or shock- reheating in the region above the 
neutrinosphere in supernovae. Hereafter, we will assume in this paper that 
is the heavy neutrino with a mass in the range 1-100 eV, and we will consider 
the two generation system of electron and tau neutrinos only. Obviously, our 
computed effects in the supernova would be identical if instead we were to choose 
Vy, as the heavy neutrino. This follows since we expect the energy spectra of 
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and Ur and their antiparticles to be nearly identical in our region of interest in 
supernovae. Working only with two neutrino generations is justified, so long as 
we assume that there exists a reasonable hierarchy of the three neutrino masses 
in which no degeneracy occurs. 

The Majorana neutrino evolution equation for two neutrino generations, in- 
cluding a vacuum mixing angle, a transition magnetic moment, and magnetic 
fields, is given by [0 











= H 

















H 



+ A sin^ 6 
|Asin2^ 




.d_ 

dr 



|Asin2^ 
+ A cos^ e 







^ + A sin^ 9 
iAsin2^ 





i A sin 26 
— + A cos^ 6 



(5) 



(6) 



where 9 is the vacuum mixing angle, and is the transition magnetic moment 
between v^. and v^-- Here is the transverse component of the magnetic field 
along the neutrino trajectory. In the usual fashion we define A = 5m? /2Ey, where 
5m? = m\—m\ > is the difference of the squared vacuum mass eigenvalues of the 
two neutrino mass eigenstates z/2 ~ v^. and vi ~ z/g, and Ey is the neutrino energy. 
We assume that the vacuum mixing angle 9 is very small, so that in vacuum with 
no magnetic fields the mass eigenstates are approximately coincident with the 
fiavor eigenstates. The effective matter potentials for z/g and Vr are given by 



ay^ = V2Gf {ue - ^Un), 
ay, = 



(7) 
(8) 



where Gp is the Fermi constant, rif. and n„ are the net number densities of elec- 
trons and neutrons, respectively. These expressions are for a neutral unpolarized 
medium and we neglect the contribution from neutrino-neutrino scattering be- 



cause its effect would be small under the conditions we consider here (see Ref. pO| 
for detailed studies of neutrino-neutrino scattering effect on MSW neutrino flavor 
transformation). In Eq. (^, Ug, rin and Bx_ are all understood to be position 
dependent in the region above the neutrinosphere which we consider here. 

By equating each two of the diagonal elements in the Hamiltonian matrix 
in Eq. (||), we flnd that there are two kinds of resonances. As expected, these 
correspond to MSW conversion and RSFP. The MSW resonance occurs when 



^Gp = A cos 29. 
The RSFP resonance occurs when 

Tie — rin) = ± A COS 29, 



(9) 



(10) 
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where the plus sign is for the Ve-^T resonance and the minus sign is for the Ve-^T 
resonance. These two RSFP resonances cannot occur at the same time. As can 
be readily seen from Eq. ([T0|), the v^-Vr resonance occurs if the sign of — Un is 
positive, whereas the t'e-^r resonance occurs if the sign of rig — nn is negative at 
the resonance position. 

The resonance region of most interest here lies above the neutrinosphere, but, 
at an early epoch (tpB ~ 0.1-0.6 s), within the radius where the shock has stalled 
(r ~ 500 km) and/or, at a later epoch (tpB ~ 3-20 s), within the radius where 
the weak and/or nuclear reactions freeze out (r ^ 40 km). Here, tpB indicates 
the time post core bounce. In the standard supernova models, Ue — Un takes 
negative values in the regions of interest. Hence, the relevant RSFP we will 
consider is ^ i^r- For a neutral medium (ng = n^), the neutron excess can be 
characterized by the electron fraction Ye, the net number of electrons per baryon. 

Ye ^ (11) 

rie + rin 

If Yf> < 0.5 then n„ > ng. From numerical supernova models, the typical val- 
ues of Ye are predicted to be about 0.4-0.45 in the region of interest above the 
neutrinosphere |2l| . 



In Fig. 1 we have plotted schematically as functions of matter density p 
the neutrino energy levels (effective mass-squared differences) in neutron-rich 
environments ( > Up) for two generations of Majorana neutrinos. To draw the 
curves in this figure we have assumed that 1/3 < Ye < 1/2. This condition on 
the electron fraction should be valid so long as we confine our considerations to 
the region well above the neutrinosphere. The number density of electrons and 
neutrons are related to matter density p in the following manner: 

ne= p YeNA, (12) 
nn= pil-Ye)NA, (13) 

where Na is Avogadro's number. The resonance density for RSFP is given by 



P^ef ^ ~ 6.6 X 10^ 



100 eV^ 



10 MeV 



cos 2^ _o 

g cm (14) 



1 - 2Ye 



whereas the MSW resonance density is 

MSW _ 1 - 21^6 RSFP f-N 
Pres — ^ Pres ■ K^'^J 
J- e 

As one can see from Eq. (plSf), the resonance density for RSFP is larger than that 
for MSW conversion as long as 1/3 < Kg < 1/2. This implies that the RSFP 
resonance takes place before the MSW resonance as neutrinos propagate from 
the neutrinosphere to the outer regions of the supernova (see Fig. 1). Since the 
matter density at the neutrinosphere is > 10^^ g cm~^, neutrinos with typical en- 
ergies and possessing cosmologically significant masses (1-100 eV) will propagate 
through one or more resonances. 
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(16) 



In order to illustrate the mechanism of RSFP, let us first work with the system 
of z/g and Vr alone, and ignore for the time being MSW flavor conversion. However, 
it should be noted that the following discussion of the effects of RSFP will be 
valid even if the MSW resonance were to coexist with RSFP, so long as the two 
resonances are well separated [see Eq. (pT]) for the non-overlapping condition for 
the two resonances]. In the limit where the vacuum mixing angle is vanishingly 
small, 6' — > 0, Eq. (||) can be reduced to the following expression, 

dr [ lyri I -f^uB a^^ + A \ [ Ur 

In this equation and hereafter we will simply write B where we have written B± 
before. However, it should be understood always that only the transverse com- 
ponent of the magnetic field along neutrino trajectories is relevant for neutrino 
spin-flavor precession. The resonance condition for RSFP in this case is given by 

V2Gf K - n,) = A, (17) 

where it is assumed that Un — > 0. At resonance, the transformation between 
Ue and Ur can be greatly enhanced even if ^uB^es <^ A, since the matter potential 
cancels the mass difference between z/g and Vr at this position [see Eq. (|l3)]. Here 
i?res is the magnitude of the transverse magnetic field at resonance. At least in a 
mathematical sense, the RSFP level crossing is very similar to the MSW one. 

The effective mixing angle 9 and the precession length L in the supernova 
environment are given by 



sm2d 



{(2/i,i?)2 + [A - V2Gf {n^ - n,)fy/^'' 
27r 



{(2/i,5)2 + [A- V2Gp (n. 
The precession length at the RSFP resonance is 



n, 



TT 



^^uBr, 



1.1 X 10^ cm 



f^i/Bxcs 



(19) 



(20) 



where fiB = e/2m is the Bohr magneton. The width of the RSFP resonance is 
given by 

(5r = 2iftan2^o, (21) 

where 



tan20n 



2.3 X 10- 



f^uBj^cs 

flB ■ G 



100 eV^ 



10 MeV 



Here H is the "density" scale height which can be expressed as 

d 



H 



dr 



ln(n^ 



-1 


d 


-1 








res 


res 



(22) 



(23) 
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To derive Eq. (p3D we have assumed that \dYe/dr\ is very small compared with 
\d\n.p/dr\. This approximation is valid in the region above the neutrinosphere 
where RSFP can significantly affect supernova dynamics and/or nucleosynthesis. 

Adiabatic spin-flavor conversion takes place when the condition Lres ^ ob- 
tains. When tan26'o ^ 1, we can employ the simple Landau-Zener approximation 
to estimate the probability for a z/^ (or v^) going through the RSFP resonance to 
remain as a z/^ (or z/g). This probability is given by 



Pi 



TT 



5r 



RSFP 



exp 



2 Lr 



exp 



-0.2 



100 eV' 
6m'^ 



10 MeV 



H 



105 



cm 



(24) 



It should be noted that the dependence of this probability on and dm? is 
opposite to the corresponding depence of the MSW survival probability [see Eq. 



It is conceivable that we may encounter situations where perhaps the mag- 
netic fields are very large, or neutrino transition magnetic moments are nearly at 
the maximum values allowed by experiment. Such a situation, in turn, could lead 
to large precession effects when tan26'o ~ 1. To properly treat the RSFP sur- 
vival probabilities in this case, we should employ the following more appropriate 
formula p2l instead of Eq. (E3), 



RSFP 



^ -exp| 



2 LroR 



cos29i cos 29 f, 



where 



cos 26',- 



A - V2GF{nn 



rip 



{(2/i,P)2 + [A - ^Gf {rin - n,)]2}i/2 



(25) 



(26) 



Here, 9i is the initial mixing angle at radius where neutrinos are produced 
(neutrinosphere), and ^/ is the final mixing angle at radius r/ where we calcu- 
late the survival probabilities. The initial mixing angle at the neutrinosphere 
always satisfies cos26'j ~ —1 for p^B{ri) < 10 /i^ ■ G, whereas cos29f ^ 1 unless 
2fii^B{rf) > A. For our choice of parameters in Sees. Ill and IV, Eqs. and 
(p5D give almost identical probabilities except around ~ 1 eV^. 

Let us now consider the case where both MSW and RSFP resonances occur 
along a neutrino trajectory. As one can see from Eq. (plSj) , the resonance densities 
for MSW conversion and RSFP will necessarily be different for Yg 7^ 1/3. This in 
turn implies that the MSW and RSFP resonances will occur at different distances 
from the neutron star. In Fig. 2 we plot a typical matter density profile at 
late times, which corresponds to tpB ~ 6 s in a numerical supernova model by 
Wilson and Mayle [^. In this figure we also indicate the resonance positions 
for RSFP (filled circles) and MSW conversion (open squares) for a neutrino with 
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El, = 25 MeV and for various labeled heavier vacuum neutrino mass eigenvalues. 
In labeling these resonance positions we have assumed that Sm"^ = rnl^ — m^ ~ m\. 
Note that the matter density at the neutrinosphere is about 10^^ g cm^^. In Fig. 
3 we plot the resonance positions (log of the radius in cm) for MSW conversion 
(dashed line) and RSFP (solid line) as functions of 5rn?. From Figs. 2 and 3, we 
can clearly see that for a given neutrino energy, the RSFP resonance occurs at 
higher density than the MSW resonace, unless > 1000 eV^. 

An analytic treatment of the case where both RSFP and MSW resonances 
occur along a neutrino path is possible if the two resonances are well separated 
in space. The non-overlapping condition for these two resonances is given by p3| 

(3Fe-l)/l"e > tan 2^ + tan 2^0- 



(27) 



If this condition is not satisfied, we would have to do a numerical integration 
of Eq. (^ to estimate reliably the survival probabilities. The above inequality 
holds for almost all of the parameter region we will discuss in Sees. Ill and IV. 
The probability for a z/g (or ) going through the MSW resonance to remain as 
a z/g (or Vj. ) is given by 



MSW 




HA sin^ 29 

10 MeV 



100 eV" 



H 



10^ cm 



sin^ 2d 
10-5 



(28) 



for 9 1. In Eq. (pSj) , we have made the approximation \d\n.ne/ dr\~^ ~ 
\d\n. p / dr\^^ ^ H. In Fig. 4 we plot the more accurate "density" scale heights 
for the RSFP (solid line) and MSW (dashed line) resonances as functions of 
rin — Tie and rig, respectively. 

By employing Eqs. (^) and (|28|), we can easily estimate the probability that 
a Vr emitted from the neutrinosphere, and subsequently propagating through an 
RSFP resonance and then later through an MSW resonance, emerges as either 
Vr, Ve-, and z/g as 



-P(^T ^t) — -Prsfp-Pmsw; 

P{^T ^e) = -pR,SFp[l — -PmSw] 
P{l^r ^Ve) = 1 - ^RSFP- 



Other probabilities such as P{ve ~^ ^r) and P(z^e 
straightforward and similar fashion. 



(29) 
(30) 
(31) 

can be estimated in a 



III. Neutrino Spin-Flavor Conversion and Hot Bubble 
r-Process Nucleosynthesis 



In the post-core-bounce evolution of the hot proto-neutron star, all six species 
of neutrinos and antineutrinos are produced thermally, provided that all vacuum 
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neutrino masses are reasonably light. These neutrinos carry away the gravita- 
tional binding energy of the neutron star on a neutrino diffusion timescale. This 
timescale is roughly ~ 10 s and is essentially set by three quantities: the mass 
of the neutron star (roughly the Chandrasekhar mass), the saturation density 
of nuclear matter, and the Fermi constant Gp- The neutrino diffusion process 
sets the timescale for all of the post-core-bounce supernova evolution. We have 
explained in the introduction that we expect the neutrinos emitted from the neu- 
trinosphere to be instrumental in heating and ejecting the envelope of material 
which surrounds the neutron star. For our purposes it is convenient, if somewhat 
artificial, to divide the evolution of this envelope into two epochs: the epoch of 
shock reheating at tpB < 1 s, and the hot bubble/r-process epoch at tpB ~ 3- 
20 s. In addition to the neutrino diffusion timescale, the neutron star radius, 
the weak and nuclear freeze-out positions, and the location of the rapid neutron 
capture environment, there is yet one more important characteristic length scale 
in the hot bubble/wind environment — the gain radius. Heating engendered by 
charged-current absorption of electron neutrinos and antineutrinos on nucleons 
wins out over neutrino losses in the hot plasma above the "gain radius," Vg > 10 
km O, m. 



As outlined in Sec. I, early deleptonization of the hot proto-neutron star 
causes its outer layers to become neutron rich. This neutron excess causes Ue to 
have a larger opacity than Ue because of the charged-current capture reactions 
on free nucleons in Eqs. and (jTB]). Consequently, have a larger average 
energy than u^. because decouple deeper in the core where the matter is hotter. 
The typical average energies for u^. and z/g at the hot bubble/r-process epoch are 
11 and 16 MeV, respectively. On the other hand, the typical average energies 
of z/^, Ufj^, and Ur are all about 25 MeV. This is because they only have the 
neutral-current reaction opacity sources common to all neutrino species. Hence, 
i/^, Ufj^, Ur and Ur decouple in regions hotter than those where z/g and z/g decouple. 
Thus, the average neutrino energies in supernovae during the post-core-bounce 
epoch always satisfy the hierarchy: 

{E^^J ^ {E,^J > {E,^) > {E,^). (32) 
The value of Y^. in the region above the neutrinosphere is determined by the 



charged-current reactions in Eqs. (p!aD and (TB). The rate of change of Ye with 
time t or radius r in this region is given by 

dY dY 

^ = vir)^ = X, - X,Ye, (33) 
dt dr 

where v{r) is the radial velocity field of the material in the supernova. In this 
equation, Ai = X^^n + Xe+n and A2 = Ai + Xp^p + Xe-p. Here, A,,^„ and Xp^p 
denote the rates of the reactions in Eqs. (p!a| ) and ([TBI), respectively, and A 



p 



and Xe+n denote the rates for their reverse reactions. At some point above the 
neutrinosphere the local material expansion rate in the hot bubble will be faster 
than the rates of the reactions in Eqs. (p!a]) and (ITH). We shall term this location 
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the weak freeze-out point, since the value of for material flow above this point 
will remain constant in time and space. Above the weak freeze-out point the 
solution of Eq. (|33|) gives 



^'^N ~ ,00 , (35) 

Evfu{Ey)dEy 



Ye{T^Y0) ~ >;(rwF0) ~ . , YBTTr 'V ^^^^ 

where rNFO and rwpo are the nuclear and weak freeze-out radius, respectively. 
Here we have neglected Ae-p and Ag+n- This is a valid approximation for our 
purposes, since the matter temperature in the region above the gain radius is 
small compared with the effective temperatures for Vf, and v^. energy distributions, 
and hence they are small compared with Ap^p and \y^n- 

The rate can be calculated as 



47rr^ 



with (z/, A^) = (z/g, n) or (z/g, p). Here, is the neutrino luminosity, jv{E^ 
is the normalized neutrino energy distribution function. We take fu{E^) to be 
Fermi-Dirac with zero chemical potential in character for all neutrino species, i.e., 

1 1 El 
^"^^"^ ^ LSOSTf l + exp[E,/r,]' ^^^^ 

where Ti, is the neutrino temperature. The temperatures characterizing the dis- 
tribution functions of each neutrino species at late epochs are approximately 
given by T^^ ^ 3.5 MeV, T^^ 5.1 MeV, and T^^^^, ^ ^ 7.9 MeV. In 

Eq. (|35D , is the cross section for the reactions in Eqs. (|l^) and ([TbD , and is 
approximately given by 

a.^^9.6x lO-^^(^)' cm^. (37) 



^ ) 



Utilizing Eqs. (|35D -(p7D, and taking into account that each neutrino species 
has roughly the same luminosity at late epochs, we can approximate Eq. (0) as 

For typical neutrino average energies (Ep^) = 16 MeV and (E^,^) = 11 MeV, 
we obtain ^ 0.41. This value can be regarded as the standard supernova 
model prediction for the Yg in the hot bubble in the absence of RSFP and MSW 
conversion. 

Let us examine now how the value of Y^ would be affected by RSFP and 
MSW resonances occurring along neutrino trajectories below the weak freeze-out 
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radius in the hot bubble. For this case, the computation of Ye must employ the 
"distorted" energy distribution functions for z/g and which will result from 
the energy- dependent flavor conversion associated with the neutrino propagating 
through an RSFP and/or an MSW resonance. By using the survival probabilities 
calculated in Eqs. (p^ and (p8|) in Sec. II at these RSFP and MSW resonances, 
we can estimate the effective and u^. energy distribution functions at the weak 
freeze-out radius to be 

/p,(E,) = /° (E,)PRSFp(i?.) + fu^Eu)[l - PRSFpiE,)], (39) 



+/0 (E,)[l - Prsfp(^.)][1 - Pmsw(^.)] 

f,^iE,)PKSFp{E.)[l - PuswiE,)]. (40) 

Here f^{Eu) represents the appropriate initial neutrino energy distribution func- 



tion. As in Eq. (|36|) , these initial neutrino energy distribution functions are all 
assumed to be Fermi-Dirac with zero chemical potential in character, but with 
different temperatures T,^^, Tp^ and Tj,^. By employing the distorted energy dis- 
tribution functions for and z/g in Eqs. (|39| ) and (|40|), respectively, we can use 
Eq. ( pSD to calculate the rate X^jy, and hence Eq. (|3^) to estimate Y^. 



We can now describe the results of our calculation of Y^ for two cases: (1) 
RSFP but no MSW conversion along a neutrino trajectory, and (2) both RSFP 
and MSW conversion on the same neutrino path. In our calculations we assume 
that the magnetic field profiles around the proto-neutron star are as follows: 

Bir) = Bo (ro/rr x lO^^ (41) 

where ro = 10 km and n = 2 or 3. A magnetic field of order 10^^ G around 
the proto-neutron star is plausible, especially given that some pulsar magnetic 
fields are at least this large. However, this argument should be viewed with some 
skepticism. There is no guarantee, for example, that the large pulsar magnetic 
fields are not generated only after the epoch of neutrino heating that we are 
investigating. Nevertheless, the RSFP effects which we describe actually depend 
on the product of the neutrino transition magnetic moment and the magnetic 
field component transverse to the neutrino trajectory. So in addition to any 
uncertainties in the magnetic field magnitude, there are geometric uncertainties 
due to the unknown distribution and orientation of the magnetic field, as well 
as the inherent uncertainties in the neutrino transition magnetic moment. For 
the purposes of our parametric study, we fix the value of the neutrino transition 
magnetic moment to be /i,^ = lO'^^/i^. Again, however, it should be noted that 
our results depend only on the combination of f^uBo for fixed value of n. 

Let us first discuss case (1) where there is only RSFP and no MSW conversion. 
This case will obtain whenever sin^ 26' <^ 10~^ and Sm"^ ~ 1-10^ eV^. We use 
the energy distributions in Eqs. (|39|) and (^0]) with Pmsw = 1 to calculate the 
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values of in this case. The initial and final mixing angles for RSFP in this 
case are taken to be at the neutrinosphere (r 10 km) and the weak freeze-out 
radius (r 40 km), respectively. In Fig. 5 we have plotted the contours of Ye (as 
labeled) in the Bq-Sitl^ plane for B oc {n = 2). These contours correspond 
to Ye values at the weak freeze-out radius, above which r-process nucleosynthesis 
may be occurring. One can see from this figure that the value of Ye is modified 
significantly from the standard no-RSFP case whenever Bo>0.1. In Fig. 6 we 
have plotted the same contours as in Fig. 7, but for B oc (n = 3). Because 
of the rapid decrease of magnetic field for n = 3, the magnitude of Bq required 
to cause similar effects on Ye is somewhat larger than that for n = 2. 

Next we examine case (2), where there exist both RSFP and MSW conversion 
along neutrino trajectories. In Fig. 7 we plot the regions in the sin^ 26-5'w? 
parameter space where neutrino flavor evolution is dominated by RSFP and/or 
MSW conversion, or neither, for the matter density proflle in Fig. 2. The solid 
and dashed lines correspond to Prsfp = and Pmsw = respectively, for 
a neutrino with Ei, = 25 MeV. We have taken fiuB^es = 1 /i^ ■ G in drawing the 
solid line in Fig. 7. This flgure clearly separates the parameter space into four 
regions, (a)-(d). In region (a), MSW conversion is nearly completely adiabatic 
and therefore very efficient in i>e ^ Vt conversion; whereas in region (b), RSFP 
is nearly completely adiabatic and thus quite efficient in Ue ^ conversion. In 
region (c), neutrino propagation is quite adiabatic through both the RSFP and 
MSW resonances; whereas in region (d), neutrino propagation through either type 
of resonance is non-adiabatic. The solid line in this figure would move upward to 
5rn? ^ 1000 eV^ if we were to take fiuB^es = 10 /i^ ■ G, and downward to 6m? ^ 3 
eV if we were to take fiuB^es = 0.1 fiB ■ G. 

In the calculations for case (2), the initial mixing angle (required for estimat- 
ing survival probabilities after neutrino propagation through an RSFP resonance) 
is evaluated at the neutrinosphere. This choice is not crucial to our results — 
it is only necessary that we evaluate the initial mixing angle in a region where 
the density is significantly larger than that at the resonance position. The final 
mixing angle ^/ is also required for estimating RSFP survival probabilities. This 
angle must be evaluated at a radius lying somewhere above the RSFP resonance, 
but below the MSW resonance. We have used the following approximation to 
estimate 9f: 



where tmsw is the radius for the MSW resonance. The non-overlapping condi- 
tion for the RSFP and MSW resonances in Eq. (|27| ) described in Sec. II must 
be satisfied if we are to employ the analytic Landau- Zener formula for survival 
probabilities. Since the two resonances are well separated so long as 6m^ < 1000 
eV^, we choose to perform calculations for this parameter range only. 

In Fig. 8 we have plotted contours of Ye (as labeled) in the sin^ 26-5rn? plane 
for the magnetic field profile of the form given in Eq. ( PD with n = 2 and 




'^MSW < '^WFO, 
'^MSW > '"WFO; 



(42) 
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Bo = (a) 0.01, (b) 0.1, (c) 1.0, and (d) 10.0 on separate plots. In plots (a) 
and (b), the magnetic field is too small to cause appreciable RSFP, and thus we 
have obtained essentially the same results as in Ref. p. However, in plots (c) 
and (d), the contours corresponding to Ye = 0.5 are significantly altered over 
the case without RSFP. In particular, in plot (d) where Bq = 10.0, there is no 
parameter region where Yf, > 0.5. We have also done a similar set of calculations 
but now with n = 3 and Bq = (a) 0.1, (b) 0.5, (c) 1.0, and (d) 10.0. Contour 
plots for these calculations are given in Fig. 9. The overall qualitative behavior 
exhibited in Fig. 9 is similar to that in Fig. 8. From both Figs. 8 and 9 we can 
conclude that the effects of RSFP will dominate over those of MSW conversion 
alone whenever Bq>1.0. Of course, this particular quantitative conclusion is 
predicated on neutrino transition magnetic moments being near their maximally 
allowed values. Obviously, smaller neutrino transition magnetic moments would 
then necessitate a larger threshold value of Bq for which RSFP would dominate 
the effects on Ye in the hot bubble. 



IV. Neutrino Spin-Flavor Conversion in the Shock Reheating Epoch 

Neutrino propagation through RSFP resonances in the region above the neu- 
trinosphere can also affect the dynamics of the supernova explosion. This comes 
about because flavor conversion in the region below the shock during the reheat- 
ing epoch, tpB ~ 0.1-0.6 s, can lead to an enhanced neutrino energy deposition 
rate in this region which, in turn, can lead to a higher shock energy. For exam- 
ple, z/g in the reheating region would become more energetic if high energy 
or Ur propagate through an RSFP resonace and become Ue- This would lead to 
an enhanced Ue capture rate on protons, and hence an enhanced heating rate. 
Therefore, in what follows, we concentrate on how RSFP and/or MSW reso- 
nances influence the charged-current neutrino and antineutrino heating rates in 
the region below the shock. 

The shock heating rate per proton or neutron is given by 

The total heating rate accompanying the Ue and z/g absorption processes in Eqs. 
(p!aD and (^) is given by 

where Yn and Yp are the number fractions of free neutrons and protons, respec- 
tively, and are approximately specifled by Eq. ( |5^ ) as F„ ~ 1 — and 1^ ~ Yg, 
respectively. The ratio of the total heating rate with RSFP (primed symbols) to 
that without RSFP (and without MSW conversion) is given by 

•/ Y'e' + Y'e'- 

^tot ^n^u„n ~l~ ^p^PeP 
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At a representative time during the reheating epoch, the temperatures of the 
relevant neutrino species are approximately given by T^^^^^ = Tp^^^^ ^ 7 MeV and 



T^^ ^ Tp^ ^ 5 MeV in the Wilson and Mayle calculations [21|. Thus, the ratio of 



the total heating rates with and without RSFP effects can be estimated to be 

T 



tot ^ Y' 



Ctot 




^ ^ 1-4, (46) 



where again we have taken the individual neutrino luminosities to be approxi- 
mately the same. This estimate of the heating enhancement factor for RSFP has 
assumed complete z/^ ^ z/g conversion in the region below the shock. Because 
Vf. and z/g have roughly the same energy distributions at the reheating epoch, 
the heating enhancement factor and the accompanying Ye value would remain 
essentially unchanged if additional MSW Vr ^ Vf, conversion were to follow the 
complete Vf, conversion by RSFP. We should note that the absolute aver- 

age energies and therefore the temperatures of the various neutrino species are a 
subject of great debate and controversy in the numerical supernova modeling com- 
munity. However, it is clear that only the differences between the temperatures 
of the relevant neutrino species are important in our estimates of the reheating 
enhancement. Neutrino transport calculation estimates of these differences are 
somewhat more reliable than those of the average neutrino energies themselves. 

In Fig. 10 we plot an example matter density profile in the region above 
the neutron star at an early epoch, tpB = 0.15 s, when shock reheating has 
commenced [?T|. At this time in the Wilson and Mayle model the shock wave is 



located at r ^ 4.7 x 10^ cm from the neutron star center. Here we employ the 
same power-law type magnetic field profiles as in Sec. HI for our estimates of the 
shock reheating enhancement, i.e., 

B{r) = Bi {ri/ry x 10^^ (47^ 

where ri = 100 km and n = 2 or 3. 

In Fig. 11 we plot the regions in the Bi-5m?' parameter space where the 
enhancement of the reheating rate is 40%. Contours in this figure are shown 
for n = 2 and 3, as represented by the solid and dashed lines, respectively, and 
are calculated from Eq. (|5]) with the use of Eqs. (||), (H), (||), and If 
the parameters (5m^ and Bi fall inside these contour lines, then the increase in 
the neutrino heating rate in the region below and near the shock is 40% more 
than that in the standard case without neutrino flavor mixing. We conclude that 
RSFP can produce heating effects which are similar to those discussed in Ref. 



25| for the MSW conversion. Therefore, RSFP may be beneficial to models of the 
supernova explosion, by virtue of increasing the average z/g energy and leading to 
a more energetic shock wave. However, we should note that the larger neutron-to- 
proton ratio (lower Ye) necessarily resulting from RSFP (as detailed in Sec. HI) 
would aggravate the problem of overproduction of the = 50 nuclei, particularly. 
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^^Sr, *^Y, and ^°Zr, which is inherent in some models of the nucleosynthesis from 
neutrino-heated supernova ejecta in this epoch 

V. Conclusions 



We have investigated the combined effects of matter-enhanced MSW conver- 
sion and RSFP of Majorana neutrinos on supernova dynamics and heavy element 
nucleosynthesis in neutrino-heated supernova ejecta. If neutrinos are of Majo- 
rana type and have transition magnetic moments of order ~ 10~^^/xb, then in 
the presence of a neutron star magnetic field of order ~ 10^^ G, resonant spin- 
flavor precession (conversion) of high energy Ur (or z/^) with vacuum masses in 
the range 1-100 eV into light Ue could increase the neutron-to-seed ratio for the 
r-process. In principle, RSFP could fix the central problem which confounds 
current models of r-process nucleosynthesis from late-time neutrino-heated su- 
pernova ejecta — obtaining a high enough neutron-to-seed ratio. We have found 
that significant enhancement in this ratio could be obtained for a plausible range 
of the proto-neutron star magnetic field, B > 10^^ G, and for neutrino transition 
magnetic moments near the maximum value allowed by the stringent astrophys- 
ical constraint. Although RSFP effects may enable r-process nucleosynthesis to 
proceed in the hot bubble, it does so at a price. The enhancement of the Ue 
energies resulting from RSFP may be at odds with the observations of these neu- 



trinos from SN1987A |13]. Of course, on the other hand, these considerations of 



the RSFP effects in supernovae may have important implications for cosmology, 
since the range of neutrino masses, 1-100 eV, required to obtain RSFP in the rel- 
evant region of the supernova is coincidently the range of interest for a significant 
neutrino dark matter component. 

We have also shown here that r-process nucleosynthesis would not be sup- 
pressed even if the MSW resonance occurs along with the RSFP. This evasion 
of the bounds from the r-process on the MSW conversion will be operative 
so long as we can be guaranteed a large value of fiuB^es > 1 /^s ■ G. For some of 
the parameters we considered, the value of Ye can be as small as 0.3 when RSFP 
for high energy neutrinos occurs near the weak freeze-out point. These very low 
values of Ye are produced in our calculations when RSFP dominates neutrino 
flavor evolution. The existence of MSW resonances is irrelevant in this case. 

We also examined the effects of RSFP during the shock reheating epoch. We 
found that RSFP can increase the total reheating rate by about 40%, but at 
the same time the concomitant reduction of Ye would exacerbate the problem of 
the overproduction of the = 50 nuclei at this epoch. This increase in total 
reheating rate may be welcome for the delayed supernova explosion mechanism, 
which relies on the energy deposited by z/g and Ue absorption reactions above the 
neutrinosphere to power the shock. On the other hand, there is at present no 
compelling necessity in supernova models for an added boost in shock energy 
from a scheme such as the RSFP of supernova neutrinos. For example, Wilson 
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and Mayle obtained a supernova explosion energy in agreement with the 
SN1987A observation by the delayed mechanism with ordinary neutrino heating 
alone. 
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Figure Captions 



Fig. 1: Schematic picture of the RSFP and MSW resonances for two generations 
of Majorana neutrinos. Each curve shows the effective neutrino mass-squared 
difference as a function of the matter density in neutron-rich (1/3 < Ye < 1/2) 
matter in the presence of a magnetic field, including effects of a neutrino vacuum 
mixing angle and a neutrino transition magnetic moment. 

Fig. 2: A typical matter density profile from numerical supernova models at late 
times (tpB = 5.8 s). Filled circles and open squares show the positions of RSFP 
and MSW resonances, respectively, for a neutrino with Ei, — 25 MeV and for the 
cases where the heavier vacuum neutrino masses are 100, 30, 10, 5, and 1 eV. 
Numbers with "eV" are for RSFP resonances, whereas numbers without "eV" 
are for MSW resonances. 

Fig. 3: Positions for RSFP (solid line) and MSW (dashed line) resonances for 
different values of Sm^. These resonance positions are for a neutrino with Ei, = 25 
MeV and correspond to the matter density profile in Fig. 2. 

Fig. 4: "Density" scale heights H as functions of n„ — (RSFP) and Ue (MSW) 
corresponding to the matter density profile in Fig. 2. 

Fig. 5: Contour plot of Yg in the Bq-5w? plane for the RSFP case where B{r) — 
Bo{ro/rf x 10^^ G (ro = 10 km). 

Fig. 6: Same as Fig. 5 but for the case where B{r) — So(ro/r)^ x 10^^ G. 

Fig. 7: Regions of the sin^ parameter space where neutrino flavor evolu- 

tion is dominated by RSFP and/or MSW conversion, or neither, for the matter 
density profile in Fig. 2. The solid and dashed lines correspond to Prsfp = 
and Pmsw = respectively, for a neutrino with E,^ — 25 MeV. We take 
IJ,uBj:es — 1 /iB ■ G in Computing the location of the solid line. 

Fig. 8: Contour plots of in the sin^ 26-5m^ plane for the cases where B{r) = 
Bo{ro/rf x 10^^ q ^ km) with Bq = (a) 0.01, (b) 0.1, (c) 1.0, and (d) 
10.0. 
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Fig. 9: Same as Fig. 8 but for the cases where B{r) = Bo{ro/r)^ x 10^^ G with 
Bo = (a) 0.1, (b) 0.5, (c) 1.0, and (d) 10.0. 

Fig. 10: A matter density profile from numerical supernova models at an early 
epoch {tpB = 0.15 s). 

Fig. 11: Regions of the Bi-Sm^ parameter space where the increase in the total 
reheating rate is 40%. The sohd and dashed hncs are for the cases where B(r) — 
Bi{ri/rY X 10^^ G (ri — 100 km) with n = 2 and 3, respectively. 



22 




23 




24 




25 




26 




Fig. 5 



27 




Fig. 6 



28 




sin^2e 



Fig. 7 



29 




30 




sin^2e sin^2e 



Fig. 9 



31 




6.6 7.0 7.4 7.8 8.2 8.6 

log ( r / cm ) 



Fig. 10 



32 




33 



